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ABSTRACT 

We present evidence from a sample of 544 galaxies from the DEEP2 Survey for evolution of the 
internal kinematics of blue galaxies with stellar masses ranging 8.0 < log M*(M0) < 10.7 over 0.2 < 
z < 1.2. DEEP2 provides galaxy spectra and Hubble imaging from which we measure emission-line 
kinematics and galaxy inclinations, respectively. Our large sample allows us to overcome scatter 
intrinsic to galaxy properties in order to examine trends in kinematics. We find that at a fixed stellar 
mass galaxies systematically decrease in disordered motions and increase in rotation velocity and 
potential well depth with time. Massive galaxies are the most well-ordered at all times examined, 
with higher rotation velocities and less disordered motions than less massive galaxies. We quantify 
disordered motions with an integrated gas velocity dispersion corrected for beam smearing (cTg). It is 
unlike the typical pressure-supported velocity dispersion measured for early type galaxies and galaxy 
bulges. Because both seeing and the width of our spectral slits comprise a significant fraction of 
the galaxy sizes, ag integrates over velocity gradients on large scales which can correspond to non- 
ordered gas kinematics. We compile measurements of galaxy kinematics from the literature over 
1.2 < z < 3.8 and do not find any trends with redshift, likely for the most part because these datasets 
are biased toward the most highly star-forming systems. In summary, over the last ~ 8 billion years 
since z = 1.2, blue galaxies evolve from disordered to ordered systems as they settle to become the 
rotation-dominated disk galaxies observed in the Universe today, with the most massive galaxies being 
the most evolved at any time. 

Subject headings: galaxies - formation, galaxies - evolution, galaxies - kinematics and dynamics, 
galaxies - fundamental properties 



1. INTRODUCTION 



In the standard picture of disk galaxy formation, 
baryons enter into dark matter halos and dissipate to 

Ttr 



form disks at their c enters (Recs fc Ostriker'1977!|White 



^Rces 1978) l^n^ Efstathiou 1980; Blumenth al et al. 



19846 . The accretion of baryons onto galaxies is expected 



to be an on-going process, and it is not yet known how 
galaxy disks respond to it, or how they evolve as mass ac- 
cretion rates change. Only recently have hydrodynamic 
simulations of galaxy formation been run with enough 
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|al."2009^'Ceverino et al. 1 20 10| [Brook et al.|201lt|Kimm et 
ial. 2011 ; Martig et al. 2012 and references therein). Such 
simulations have yet to address the evolution of internal 
galaxy kinematics over a significant period of time, which 
is the cleanest observational probe of disk galaxy assem- 
bly. Internal kinematics tell us directly about the dynam- 
ical state of galaxies, reveal the potential well depths of 
individual galaxy-dark matter halo systems, and can be 
measured for a large sample of galaxies over a significant 
range in redshift. 

We have acquired a wealth of information about blue 
galaxies over the last ^ 8 billion years since z — 1.2 
which can be roughly divided into two categories: those 
which imply a significant amount of evolution and those 
which do not. Among the observations which do not im- 
ply much evolution are those of galaxy luminosity func- 
tions. Luminosity function studies find that the num- 
ber density of blue galaxies is unchanging, and that blue 
galax ies fade by only ^ 1 i?-band magnitude since z ~ 1 



Willmcr et al. 2006; Faber et al. 2007 Bell et al 



2004 1 . Similarly, studies of stellar mass functions of blue 



galaxies find no evolution to wi t hin uncertaintie s (Bundy 



let al. 



20061 



I Borch e t al.1'2006 

though see |Drory~ al. 200!^ 
are only marginally smaller at z 



06) |P( 



P ozzetti et al.|)2dlO| al 

addition, galaxy sizes 
= 1 compared to local 
galaxies (by a factor of 1.4; Button et al. 2011). How- 
ever, despite all this evidence for a slowly evolving pop- 
ulation, there are some strong indications of substantial 
evolution. Compared with blue galaxies of similar stellar 
masses today, those at z ~ 1 have significa ntly higher 
star- formation rates by a factor of ~ 10 (e.g., Noeske et 
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Kassin et al. 



al. 20071, higher molecular gas fractions (for at least a 
tew ot the brighter galaxies which have been observed so 
far ) by a factor of 2 — 4 ( Daddi et al.||2010] Tacconi et 
al. 2010|), and inore disturbed morphologies (e7g., [Abra- 
ham et al. 1996 Abraham & van den B ergh ioOl]. In' 



summary, blue galaxies do not evolve much since z ~ 1 
in luminosity, stellar mass, or size, but they do evolve 
strongly in star- format ion rate, (likely) molecular gas 
fraction, and morphology. In this paper, we attempt to 
understand this discrepancy and the mechanisms behind 
it by studying the internal kinematics of blue galaxies 
since z ^ 1. 

Due to the large intr insic scatter in the kinematic prop- 
erties of blue galaxies ( Weiner et al.|2006a|b Kassin et al. 

hich is r 



20071, a sizable samp' 



(> 100) wh 



representative 



m terms of galaxy properties is required to study inter- 
nal galaxy kinematics over a significant look-back time. 
The first large study to address the kinematic evolution 
of blue galaxies since z ~ 1 was Weiner et al. (2006a b| 
which used data for 1089 galaxies from the i'KRS Survey. 
Among other things, they found that about one-third 
to one-half of emission-line galaxies over 0.1 < z < 1.5 
have a significant or dominant component of disordered 
motions, as measured via an integ rated gas velo c ity dis 



persion {ag] as discussed in §4.1). Weiner et al. (2006a I 



proposed a new velocity indicator to trace galaxy poten 
tial well depths. It in corporates both rotation velocity 

Kassin et al. (2007) 



(hereafter K07) followed up on 

Weiner et al._ (2006a b). Their dataset consisted of 544 
galaxies from the UEEP2 Survey over 0.1 < z < 1.2 with 
medium-resolution spectroscopy from which kinematics 
were measured and Hubble/ ACS imaging. The Hubble 
images were used to measure galaxy morphologies and 
axial ratios (which were used by K07 to correct rotation 
velocities for inclinations of galaxies to the line-of-sight). 
The main results of K07 are: 

• K07 showed that scatter in the stellar mass Tully- 
Fisher relation (TFR), i.e. the relation between 
galaxy stellar mass and Vrot, is mostly to low Vrot 
and is a strong function of morphology. Galaxies 
which scatter to low Vrot generally have disturbed 
or compact Hubble morphologies. Similarly, from 
samples of 33 an d 68 galaxies i n the IMA GES Sur- 
vey at z ~ 0.6, Flores et al. (2006) and Yang et 
al. (2008|), respectively, found a large scatter in the 
11 R due to galaxies with "perturbed or complex 
kinematics." Most previous studies of Tully-Fisher 
at both low and high redshift tried to minimize 
scatter by selecting galaxies whi ch are morpho 
logically well-ordered disks (e.g., fVerheijen 1997 



Kassin, de Jong, fc Weiner 2d06[ ~ PTzagno et al. 
20071 and references therein for low redshift, and 



e.g., 'Consclic e et al.||2005l [Miller et al.1|20TT] and 

references therein for z <; l.Sj! 

• K07 demonstrated that galaxies which scatter to 
low Vrot in the TFR have a significant contribu- 
tion to their kinematics from disordered motions 
(as quantified by <7g). 

• K07 showed that when Vrot and ag are combined 
into a velocity indicator created to trace galaxy po- 
tential well depths {So,5 as defined above), a tight 



relation with stellar mass results. This relation is 
independent of galaxy morphology, non-evolving 
since z = 1.2, and coincident with the Faber- 
Jackson relation for early-type galaxies. It demon- 
strates that galaxies are likely approximately viri- 
alized since z ~ 1 in that they have all the en- 
ergy they will have today already at z ~ 1. This 
tightening of Tully-Fisher once S0.5 is adopt ed was 
later found by other observational studies ([Cresci 
et al. "2009'; Puech et al. '2010'; 'Lem oine-tjusseroIIe 
~h Lamarcillc 2010; Lc moinc-Busscro lle et al.[ [2010[ 
"lycrgani ct al. 2012; C atinella et al.j 2012) a ndnu^ 
merical simulations of interacting galaxies ( ,Coving- 
ton et al.[[20T0| . 



• K07 found that the majority of scatter in the TFR 
at low redshift is at lower stellar masses (< 10^° 
M0). They found that the scatter at higher masses 
was low at low redshift, and increased with increas- 
ing redshift to z = 1.2. K07 hypothesized that 
this meant higher mass galaxies have been settling 
onto the TFR since z = 1.2, and that lower mass 
galaxies are settling onto the TFR today. Simi- 
larly, studies with the IMAGES Survey argue for 
significant kinematic evolution since z ^ 0.6 (e.g., 
Yang et al.||2008l [Neichel et aL][2008l [Puech et al.| 



2008). 



The K07 study leaves a number of open issues: 



• First, although it is clear that galaxies are evolving 
in Vrot anc 
quantified. 



m Vrot and ag since z = 1.2, this evolution is not 



• Secondly, it is also clear from K07 that kinematic 
evolution is a function of galaxy stellar mass, but 
this is not quantified or addressed in detail. 

• Finally, although K07 found that the relation be- 
tween 5*0.5 and stellar mass relation does not evolve 
to within uncertainties since z = 1.2, it has yet to 
be determined whether galaxies evolve in 5*0.5. 

In this paper we address these open issues with a more 
sophisticated analysis of the dataset in K07. In §2 we 
detail how the K07 sample was selected; this was not 
included in K07 because it was a letter paper and did 
not have the space to do so. In §3 we discuss the stellar 
mass and Hubble size measurements which are adopted 
from other studies, and describe the measurement of 
emission line extents made in this paper. In §4, since 
it was not addressed in detail in K07, we describe how 
Vrot and CTg are measured from emission lines in galaxy 
spectra and discuss what they measure. We then show 
how V^oijCg) and 5o.5 evolve with redshift in §5, and 
discuss the influence of galaxy stellar mass and size on 
this evolution in §6. In §7 we quantitatively define a 
settled disk galaxy and study how the fraction of settled 
disk galaxies evolves with redshift. Our findings are com- 
pared with previous measurements of galaxy kinematics 
in the literature over < z < 3.8 in §8. Conclusions are 
given in §9. A ACDM cosmology is adopted throughout 
{h = 0.7, 57„i — 0.3, r^A = 0.7), and all logarithms are 
base 10. 
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Fig. 1. — Rest-frame U — B color (AB system) versus stellar 
mass diagrams for galaxies with solid spectroscopic redshifts in 
field 1 of the DEEP2 Survey (grey points) and for the K07 sub- 
sample which is used in this paper (black points) are shown in bins 
of redshift. The K07 sample provides good coverage over galaxy 
mass. Although the K07 sample is moderately biased towards bluer 
galaxies at a fixed mass, the effect is not severe. Dotted lines 
demarcate the mass-limited sample used later in this paper. 

2. SAMPLE SELECTION 

In this paper we use the dataset of K07. Since K07 
was a letter paper, it did not include a detailed descrip- 
tion of how its galaxy sample was selected, so we include 
it here. One of the main differences between K07 and 
most other studies of galaxy kinematics lies in its sam- 
ple selection. Typically, constraints on galaxy morphol- 
ogy are applied at the outset to remove disturbed and 
major-merger galaxies. K07 instead included all galaxies 
with emission lines which are bright enough to measure 



kinematics from, regardless of morphology. Since dis- 
turbed morphology and disordered motions (as measured 
by ag) are correlated, removing galaxies with disturbed 
morphologies would introduce a strong bias toward only 
the most rotation-dominated systems (K07). Little or no 
evolution in internal galaxy kinematics would be found 
for a sample selected in this manner. Moreover, such a 
sample would not be representative of galaxies at higher 
redshifts s ince they tend to have disturbed morpholo 
gies (e.g ., iGrifhths ct al.||1994' lwT7Ti3:7::rT:r-TTrT7T7TF 



Abraham et al., 1996 



Windhorst et al."1995 



1 



^ Abraham &z van den Bergh 2001 

Even when a strIct~morphological selection is loosened 
somewhat to include galaxies which are moderately dis- 
turbed but disk-like, still little scatter and no evolution 
in Tully-Fi sher since z = 1.7 are found (Miller et al. 
2011 2012). Furthermore, many studies of galaxy kine- 
matics at medium to high redshift are limited to galaxies 
with extremely bright emission lines, which are not rep- 
resentative of the general population of galaxies. The 
sensitivity of DEEP2 allows the K07 sample to probe a 
large range in galaxy stellar mass and color (Figure [T]) . 

In addition to broad morphological selection criteria 
and a large range in stellar mass, in order to study the 
evolution of disk galaxy kinematics over a significant 
look-back time, medium resolution spectra of a sizable 
sample of galaxies (> 100) which span a wide range in 
redshift are required. A sizable sample is needed because 
emission-line galaxies have significant intrinsic sca tter in 

|2006a|b 
large co 



their kinematic properties (Weiner et al 
The K07 sample derives trom 



K07). 

lection 



a very 

of 9715 galaxies with solid spectrosc opic redshifts from 



Keck/DEIMOS ( Faber et al.| |2003[ ) in field 1 of the 
DEEP2 Survey (May 20U7 catalog), as described below. 
Field 1 is home to the AEGIS multi-wavelength survey 



pavis et al .||2007|. The hmiting magnitude of DEEP2 
is Hab < 24.1 and galaxies in field 1 are not subject to a 
color selection ( ,WiUrner et al.|[2"006 Newman et al. sub- 
mitted) . In Figure |l] we show rest-trame color-magnitude 
diagrams for galaxies in field 1 of DEEP 2 and highlight 
the K07 sample. The K07 sample provides good cover- 
age over galaxy mass. Although it is moderately biased 
toward blue galaxies at a fixed mass, the effect is not 
severe. 

In the following we detail how galaxies were selected 
from field 1 of the DEEP2 Survey to be in the K07 
sample. First, although DEEP 2 reaches a redshift of 
z — 1.5, K07 chose to be conservative and limit their 
sample to z = 1.2, leaving 8787 galaxies. This helps 
to avoid the high-redshift end of the survey where the 
selection is in the rest-ultrayiolet and becomes depen- 
dent on color ( Willmer et al.|[2006 |. It also ensures that 
galaxy axis ratios, which are used to correct rotation ve- 
locities for inclinations of galaxies to the line-of-sight, are 
measured from Hubble imaging which has not yet shifted 
into the ultraviolet since these wavelengths may not trace 
the bulk of the stellar mass of galaxies. Next, galaxies 
were chosen to have Hubble Space Telescope/ ACS imag- 
ing at V and /, which further reduces the sample to 
3523 galaxies. This is because there are portions of field 
1 which are not imaged by Hubble. Aside from removing 
galaxies which did not have Hubble imaging, this require- 
ment did not further affect our sample selection because 
the Hubble images are deeper than the spectral survey. 
For our most important selection, galaxies with bright 
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Fig. 2. — Relationshi ps between kinem atic quantities (erg, Vrot, So. 5) and stellar mass are shown in bins of redshift for the K07 sample. 
Stellar masses are from |Lm et al.|2007| The Vrot and So. 5 relations were shown previously in K07, and we show the ag relation here for 
completeness. While large scatter is tound for the ag and Vrot relations, the So. 5 relation is significantly tighter and does not evolve to 
within uncertainties since 2 = 1.2. The fit to the So. 5 relation for the lowest redshift bin is shown as a solid line (slope of 0.30it0.2, intercept 
of 1.93 it 0.01 at log Af, = 10 Mq ; K07), and is repeated in the higher redshift bins as a dashed line. It is consistent with the Faber- Jackson 
Relation for early type galaxies, where a is measured from stellar absorption lines (K07). Hashed regions denote areas where the survey 
is not sensitive, and dotted lines demarcate the mass-limited sample used later in this paper. If a galaxy is deemed too morphologically 
disturbed to accurately determine an inclination (9% of the sample), Vrot is not inclination corrected and the galaxy is shown as an upper 
limit symbol in the Vrot and 5o.5 plots and an open circle in the ag plot. 



enough emission lines to measure kinematics (> 10 
erg cm^) were chosen, leaving 1692 galaxies. 

Galaxies were further selected by K07 to have inclina- 
tions measured from Hubble images between 30° < z < 
70° , and spectrographic slits with position angles aligned 
to within 40° of their major axes. These cuts leave a 
sample of 755 galaxies. Inclinations are measured from 
V + /-band H ubble images using the SExtractor soft- 
ware package ([Bertin fc Arnouts 1996), as described in 
Lotz et al. ( 200(jp , and have an uncertamty of ^ 10°. The 
inclination requirement avoids nearly face-on galaxies for 
which inclination measurements are very uncertain, and 
minimizes the effects of dust in the determination of stel- 
lar mass for highly inclined systems. However, if the mor- 
phology of a galaxy was deemed by eye to be disturbed 
enough such that its inclination could not be reliably de- 
termined, it was included in the sample but not corrected 
for inclination. There are only 24 of these galaxies and 
they are flagged in most of the following analysis. Re- 
moving them from the sample does not significantly af- 
fect the results of this paper. For the slit position angles, 
due to the large slit width (1") compared to the apparent 
sizes of galaxies in our sample 3") and the effects of 



seeing, accurate kinematics can be measured if slit po- 
sition angles are of fset by up to 40° fro m galaxy major 
axes (Figure 13 of Weiner et al. 2006a). Finally, if an 
emission line was affected by a sky line or an instrumen- 
tal artifact, the corresponding galaxy was removed from 
the sample, leaving a total of 544 galaxies in the K07 
sample. No constraints on morphology were applied. 

To study the evolution of galaxy kinematics for a sam- 
ple of blue galaxies which is mass-limited at all redshifts, 
in §5.1 we limit the K07 sample to galaxies with stellar 
masses over 9.8 < log Mh,(M0) < 10.7, as demarcated 
in Figures [T]~ [3) Since the stellar mass function of blue 
galaxies does not evolve significantly since z ^ 1 (|Borch 
et al.||2006[ |Bundy et ah, 2006; Pozzetti et al ||2010p , we 
do not change this stellar mass range with redshift. This 
cut removes 274 galaxies from the K07 sample, leaving 
270. Throughout the paper we will refer to this as the 
"mass-limited sample." 

3. STELLAR MASS AND SIZE MEASUREMENTS 

In this paper we will study the relationship between 
galaxy kinematics and stellar mass, building upon the 
scaling relations st udied in K07. W e adopt stellar mass 
measurements from Lin et al. ( 2007 ) via private commu- 
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Fig. 3. — Left: The ratio of ordered to disordered motions in a galaxy, Vrot/o-g, is shown as a function of stellar mass and in bins of 
redshift for the K07 sample. Galaxy symbols are as in Figure[2] Hashed regions denote areas where the survey is not sensitive, and dotted 
lines demarcate the mass-limited sample used later in this paper. The galaxies define an upper envelope which does not evolve with redshift. 
However, galaxies evolve within the envelope: the lower-right corner evacuates with time (from high to low redshift), with higher mass 
galaxies evolving the fastest. Right: This evolution is more clearly shown as histograms in Vrot/fg for galaxies divided into stellar mass 
bins. 

nication, as was done in K07. The uncertainties on these 
stellar mass measurements are ~ 0.25 dex. 

In addition to the effects of stellar mass on galaxy kine- 
matics, we will also look at galaxy sizes. Two size mea- 
surements are studied: emission line extent and contin- 
uum size. In this paper we measure emission line ex- 
tents from ground-based DEEP2 spectra from the same 
lines used to measure kinematics. However, unlike the 
kinematic measurements, no seeing correction is applied. 



The spatial extent of the line emission is measured as 
foll ows for a given galax y. For details on the method, 



see 



Weiner et al. ( 2006a I . First we measure the spectral 



continuum by collapsing the two-dimensional spectrum 
in wavelength over a 15A range around the line and two 
lOOA spatial ranges on either side of the line. We then fit 
a Gaussian profile along the slit using a non-linear least- 
squares routine and adopt the resulting full-width-half- 
maximum (FWHM) as our measurement of emission line 
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Fig. 4. — For the mass-limited sample, the distribution of galaxy 
stellar masses among the redshift bins in Figure |3] is shown. It 
demonstrates that the average stellar mass of the sample is not the 
same at all redshifts: it is lower at lower redshift. This is likely 
because higher mass galaxies are leaving the sample to become 
early types. This phenomenon works in the opposite direction to 
our results, making them lower limits to the intrinsic evolution, as 
explained in §5.2. 

extent. For the galaxy continuum sizes, we adopt mea- 
surements made f rom Hubble images usi ng the GIM2D 
softw are package (Simard et al. 2002) by Button et al. 
|2011| and E. Cheung et al. in preparation, 'i'o take into 
account band-shifting with redshift, sizes are measured 
from V and /-band images for galaxies at redshifts less 
than and greater than z = 0.6, respectively. Since the 
continuum sizes are measured from Hubble images, they 
are unaffected by seeing. 

4. KINEMATIC MEASUREMENTS 

Measurements of internal galaxy kinematics from emis- 
sion lines are adopted from K07. In this section we dis- 
cuss some relevant details and the phenomena they de- 
scribe. Further det ails on how these measu rements were 
made are given in Weiner et al. ( 2006a|b ) for a similar 
survey with a lower spectral resolution; R ~ 2100 com- 
pared with R ^ 5000 for the K07 sample. 

For each galaxy the emission line used to measure kine- 
matics is the one with the highest signal-to-noise in its 
spectrum. For the vast majority of galaxies, the Ha 
A6563, [O II] AA3726.0, 3728.8, and [O III] A5007 lines 
are used. Rotation velocities on the flat parts of the ro- 
tation curves (Vrot x sm(i)), where i is the inclination 
of a galaxy to the line-of-sight) and integrated gas veloc- 
ity dispersions (tTg) were measured simultaneously from 
the emission lines, and the effects of seeing (typically 
^ 0.7") were taken into account. Our spectral resolution 
allows for Vrot x sin{i) and (Tg to be measured down to 
^ 5 and ~ 15 km s""'^, respectively. An uncertainty of 
10 km s~^ is adopted for both Vrot x sin{i) and Ug to 
account for random errors and the dependence of model 
parameters on the assumed seeing and scale radius of the 
rotation curve. Although many of the rotation curves do 
not "turn over" due to seeing, our model is still able to 



fit for the r otation velocity on t he flat part of the rota- 
tion curve (Weiner et al. 2006a §2.3.2). In addition, it 



has been demonstrated using deeper data with the same 
telescope/instrument which observe a turn-over in the 
rotation curve for 90% of galaxies that spect ra of our 
depth do not show a bias in Vrot (Figure 6 of IMiller et 
Benefits of the deeper data (as used in Miller et 



al. 



2011) 



2011 ) are smaller errors in Vrot x sin(i) (and therefore 



less scatter) and the ability to probe galaxies with fainter 
emission lines. Except for the 24 galaxies for which reli- 
able inclinations could not be determined (§2), values of 
Vrot X sin{i) are corrected for galaxy inclinations mea- 
sured from the Hubble images. The results of this paper 
do not change significantly if these 24 galaxies are re- 
moved from the sample. 



4.1. What Gg and Vrot Measure 

The integrated gas velocity dispersion we measure (<Tg) 
is unlike the typical pressure- supported velocity dispersion 
measured for early type galaxies or galaxy bulges. This is 
the case for two reasons. First, ag is measured from emis- 
sion lines which trace the gas in galaxies, as opposed to 
absorption lines which trace stars. Gas, unlike the colli- 
sionless stars in an early type galaxy, can dissipate energy 
and therefore cannot remain in a high dispersion equilib- 
rium state with crossing orbits. Secondly, because the 
spectral slits used in the DEEP2 Survey are wide com- 
pared to the apparent sizes of the galaxies observed (1" 
slit vs. ~ 3" galaxies), (jg effectively integrates velocity 
gradients on scales at and below the seeing limit (|Weiner| 
et al.||2006a] [Covington et al.||2010| [Epinat et aT.pOlOp . 

35 km s~^, which is the upper 
limit for well-ordered disk galaxi es in the local Universe 



(Figure 15 of|Epinat et al. 2010), ag measures the rela- 
tive motions oi individual star-torming regions in galaxy 
disks. However, for values of Ug greater than ^ 35 km 
s~^ , (Tg integrates over velocity gradients which can corre- 
spond to non-ordered gas kinematics such as small-scale 
velocity gradients, gas motions due to star-formation, or 
superimposed clumps along the line of sight. In this pa- 
per, we collectively refer to these velocity gradients as 
"disordered motions. " 

To test for interference from Vrot in measurements of 
(Tg, in §5 we compare Og measured for galaxies in our 
mass-limited sample with galaxies from DEEP2 in the 
same stellar mass range which are face-on in Hubble im- 
ages (i.e., i < 30°). We find no significant difference in 
the trends of ag with redshift between these two samples, 
making it unlikely that the rotation velocity of a galaxy 
interferes with a measurement of its ag. 
To further inyestigat e what ag and Vrot measure, in 



Covington et al. (2010) we performed mock observations 
ot a suite ot simulations of major mergers of disk galax- 
ies. These simulations provided model galaxies in various 
stages of disorder, ranging the gamut from well-ordered 
disks to merger remnants. The mock observations were 
performed just as the actual observations. Galaxies were 
redshifted to z ~ 0.3 and z ~ 1.0 and seeing, slit width, 
and detector pixel size were taken into account. Our al- 
gorithm for fitting kinematics was found by |Covington| 



et al.| (2010) to be successful at reproducing ag, but was 
tound to underestimate Vrot by up to 30%, independent 
of merger stage or redshift. However, as we note above. 
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Miller et al. (20111 did not find a systematic difference 
in Vrot between the deeper data and data similar to ours. 
Therefore, if Vrot is underestimated by up to 30%, it is 
probably underestim ated by this much in both our sam- 
ple and in the deeper Miller et al. (2011) dataset. 



5. EVOLUTION OF GALAXY KINEMATICS 

In this section we examine how the kinematics of blue 
galaxies evolve with redshift. First we examine the full 
K07 sample, and then wc study a mass-limited sub- 
sample. We return to the full sample again in §7. 

5.1. The Full K07 Sample 

In Figure [2] we show Cg, Vrot, and 5*0.5 versus stellar 
mass in bins of redshift for the K07 sample. The Vrot 
and 5*0.5 relations were already shown in K07; here we 
show the (Tg relation for completeness. The Vrot versus 
stellar mass relation (i.e., the stellar mass Tully-Fisher 
Relation) shows large scatter to low Vrot at all redshifts 
(K07), as does tXg versus stellar mass. However, when 
Vrot and a a are combined into 5 *0.5 (a tracer of potential 



well depth; Weiner et al 



,|2006a[ 

mtens s 



K07) , the resulting rela- 
tion with stellar mass tightens significantly and does not 
evolve with redshift since z = 1.2 to within uncertainties 
(K07). Figure I2] also shows that at lower redshift higher 
mass galaxies nave on average larger values of Vrot and 
smaller values of a-g than at higher redshift (K07). It 
is also clear from this figure that at lower redshift lower 
mass galaxies have low values of Vrot and high values of 
(Tg (K07); they are too faint to be observed at higher 
redshifts in our sample. 

These trends are illustrated more clearly in Figure 
(left column) where the quantity K-ot/fg, the ratio 
ordered to disordered motions in a galaxy, is plotted ver- 
sus stellar mass in bins of redshift The galaxies define 
an upper envelope in Vrot/o'g versus stellar mass-space 
which does not evolve since z = 1.2. However, galaxies 
evolve within the envelope: with increasing time (i.e., 
decreasing redshift) they move away from the bottom of 
the plot toward the upper envelope, and this happens 
first for massive galaxies. This evolution is quantified by 
the histograms in Figure |3] (right column) which demon- 
strate that galaxies move towards higher Vrot I with de- 
creasing redshift (i.e., with time), and the highest-mass 
galaxies have on average the highest values of Vrot / Cg at 
all times. 

5.2. A Mass-Limited Galaxy Sample 

To investigate these trend of Vrot/<^g with redshift for 
all but the very lowest and very highest mass galaxies in 
our sample, we use a mass-limited sample of galaxies with 
stellar masses over 9.8 < log M*(M0) < 10.7 (vertical 
dotted lines in Figures l]|3l, as described in §2. For the 
remainder of this section and for §6, we will focus on this 
mass- limited sample of 270 galaxies. In Figure |4) for the 
mass-limited sample we plot the distribution of galaxy 
stellar masses among the redshift bins in Figure [3] It 
shows that the average stellar mass is not the same for 
all redshift bins: at lower redshifts it is shifted towards 
lower masses. This is likely because higher mass galaxies 

In the evaluation of Wot/fj, a minimum value of Ug of 15 km 
, which is the lowest value which we can measure, is adopted to 
avoid dividing by small numbers with relatively large uncertainties. 



are leaving the sa mple as they tran sform to red early 
type galaxies (e.g., Faber et "al]|2007 ) . This trend works 
in the opposite direction to our results, as we will discuss 
later in this section. Therefore, the evolution we find is 
interpreted as a lower limit to the intrinsic evolution. 

Figure [51 shows the evolution of ag, Vrot and 5o.5 with 
redshift for the mass-limited sample. All three quantities 
have large intrinsic scatter. We will demonstrate that a 
portion of the scatter is due to stellar mass or size in §6. 
Medians of the individual measurements in bins of red- 
shift are shown as red triangles in Figure [5] Errors on 
these medians are calculated by bootstrap re-sampling 
the data in each redshift bin. Bootstrap resampling is 
used as a standard method for estimating the statistical 
error without assumptions about the shape of the distri- 
bution. We re-sampled with replacement from the data 
1000 times, and calculated the dispersion of the sample 
medians. This dispersion is adopted as the error on the 
median points. It accounts for the spread in the indi- 
vidual data points due to both intrinsic variation and 
observational error. The errors on the medians are small 
because the sample is large and there are many points in 
each redshift bin. 

We fit linear relations to ag, Vrot and 5*0.5 versus red- 
shift by performing linear least-squares fits to the binned 
medians, taking into account the errors on the medi- 
ans. The results are not sensitive to the redshift bins 
used. We fit to the medians rather than directly to the 
data because the distributions of the data points are non- 
Gaussian and the intrinsic spread of the points is signifi- 
cant compared to the individual error bars, which would 
cause standard least squares fits to be misleading. To 
avoid covariances on the fitted parameters, the medians 
are zero-pointed near the middle of the sample when per- 
forming the fits. Of the three quantities, the median ag 
shows the strongest evolution (5. Oct). It increases linearly 
with redshift to z = 1.2 (i.e., decreases with time) as 

logCTg-0.42= (0.25±0.05)(z-0.82) + (1.16±0.01), (1) 

and the fit has a od 1.1. The median ag increases from 
27 ± 1 km s^i at z = 0.2 to 47 ± 1 km s^^ at z = 1.2. 
The quantity which shows the second-strongest evolution 
is the median Vrot (4.2ct). It decreases with redshift as 

log Vrot -2.00 = (-0.21 ±0.05)(z- 0.82) -(0.04± 0.02), 

(2) 

and the fits has a of 16.1. The median Vrot decreases 
from 123±1 km s'^ at z = 0.2 to 76±1 km s^^ at z = 1.2. 
Finally, the median 5o.5 shows the weakest evolution of 
the three kinematic quantities (3.6ct). It decreases with 
redshift as 

log 5o5-l-92 = (-0.11±0.04)(z-0.82)-(-(0.00±0.01), 

(3) 

and the fit has a of 2.8. The median 5o.5 decreases 
from 97±1 km s'^ at z = 0.2 to 75±1 km s'^ at z = 1.2. 

Some of the evolution found may be due to a changing 
distribution of galaxy masses in the sample with time. 
This can be caused by galaxies growing in mass and/or 
high mass galaxies quenching their star-formation and 
leaving the sample. To address this we remove the de- 
pendence on potential well depth by examining ratios 
with 5*0.5, namely Vrot/So.5 and ag/So,^. These ratios 
are shown as a function of redshift in Figure [6j They 




Fig. 5. — The evolution of ag, Vrot, and So. 5 = \JQ.5V^^^ + (Tg with redshift for the mass-limited sample (9.8 < logM,(MQ) < 10.7) are 
shown. Galaxies increase in ag and decrease in Vrot and 5o.5 with increasing redshift. Galaxies are plotted as black symbols, as described 
in Figure |2] Binned medians are shown as red triangles with error bars representing the error on the medians determined by bootstrap 
re-sampling the data. Linear fits to the binned medians are shown as red lines. Medians errors on the individual measurements are shown 
in the lower right-hand corners of the plots, and include uncertainties in galaxy inclinations for Vrot and 5o.5. 
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Fig. 6. — The evolution of Ug/So.s, Vrot/ So, 5, and VVot/fg with redshift for the mass-limited sample are shown. The first two of these 
ratios remove the dependence on potential well depth (i.e., So. 5) in the ag and Vrot plots in Figurejs] Symbols, error bars, and fits are the 
same as in Figure |5] 



express the fraction of kinematic support which comes 
from ordered and disordered motions. The median of 
the ratio Vrot/5'o.5 decreases with redshift (3. Oct) as 

log K-ot/^o 5-0.10 ^ (-0.06±0.02)(z-0.82) + (0.00±0.01), 

(4) 

and the fit has a of 4.5. The median of the ratio 
<Tg/Sa,^ increases with redshift (5. Oct) as 

log (Tq/Sa 5+0.35 = (0.45±0.09)(z-0.82) + (0.01±0.03), 

(5) 

and the fit has a of 2.3. As before, these are lin- 
ear least-squares fits to the median points. These ratios 
demonstrate that the increasing role of ag and the declin- 
ing role of Vrot with increasing redshift are not due to a 
changing distribution of galaxy masses in the sample. 

The evolution found is interpreted as a lower limit to 
the intrinsic evolution. As discussed earlier in this sec- 
tion and shown in Figure |4j the average stellar mass 
of our mass-limited sample is lower at lower redshifts. 
From Figure |3] it is apparent that on average low mass 
galaxies have higher ag and lower Vrot than more massive 
galaxies at any redshift (we will investigate this further 
in §6). Therefore, our findings that for galaxies in the 



mass-limited sample the average ag and Vrot are lower 
and higher, respectively, at lower redshift than at higher 
redshift are interpreted as lower limits to the intrinsic 
evolution. 

In addition to ratios with 5*0.5, we look at Vrot/cg, a 
ratio of ordered to disordered motions^'^. This is not an 
ideal quantity since neither Vrot nor ag is constant with 
redshift. Nevertheless, due to the significant trends of 



ag and Vrot with redshift, the median Vrc 
strong evolution (4.8(7; Figure [6|: 



also shows 



log Kot/cTg- 1.55 = (-0.48±0.10)(z-0.82)-(1.13±0.04) 

(6) 

with = 2.8. 

To determine whether the Vrot of a galaxy influences 
the measurement of ag, in Figure [Tj we examine the evo- 
lution of ag with redshift for galaxies which are face-on 
in Hubble images (i.e., i < 30°) and in the same stel- 
lar mass range as the mass-limited sample. This is a 
subset of galaxies for which ag should be completely se- 
cure. We find the same evolution in the median ag to 
within uncertainties for face-on galaxies: log ag — 0.42 = 
(0.24±0.04)(z-0.82)-h(1.20±0.02), demonstrating that 
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Fig. 7. — The evolution of ag with redshift for galaxies oriented 
face-on in Hubble images (i < 30°) (black points, red triangles show 
medians in redshift bins, solid red line shows fit to median points) 
is compared with that for galaxies in the mass-limited sample from 
Figure [5] (dashed red line shows fit). Face-on galaxies provide a 
pure measurement of tTg, uncontaminated by Vrot- Both samples 
span the same stellar mass range, and both fits are given in the text. 
There is no significant difference between the fits, demonstrating 
that our measurements of ag are unaffected by galaxy rotation. 

Vrnt does not influence the measurement of a„. 



mass bin differ significantly from the other mass bins. 

Next we examine trends with size. We consider two 
size measurements: Hubble continuum sizes and emis- 
sion line extents, as described in §3. Hubble sizes are 
unaffected by seeing, and emission line extents are af- 
fected by seeing and are not corrected for it. Figure |9] 
is similar to Figure |8] except galaxies are color-coded ac- 
cording to size instead of stellar mass. For both size 
measurements, three size bins are shown: < 4 kpc, 4-6 
kpc, and > 6 kpc. Figure |9] shows that on average at 
all redshifts examined the largest galaxies are the most 
kinematically well-ordered, the smallest the most disor- 
dered, and intermediate-sized galaxies have properties in 
between the two. This is the case for both size mea- 
surements. It demonstrates that emission lines, which 
we use to measure galaxy kinematics, are likely good 
tracers of the spatial distribution of the galaxy continua 
observed in Hubble images. It also demonstrates that we 
are able to measure reliable rotation velocities for even 
the small galaxies in our sample. Similar to the stellar 
mass bins, there is also a size threshold. It is at 4.0 kpc, 
which corresponds to the threshold in stellar mass (log 
— 10.4 M0) in the size versus stellar mass relation 
for the mass-limited sample. 

Figures |8] and |9] also demonstrate that the evolution in 
Figures [5] and [6] is not due to a special population of high- 
redshift galaxies with high values of ag and low values 
of Vrot- It occurs even if only massive/large galaxies are 
considered. 



6. DOWNSIZING 

In this section we demonstrate that much of the scat- 
ter in the relations of (Jg/Sd^ and Vrot/So.5 with redshift 
is attributable to galaxy properties such as stellar mass 
and size, and find that trends with stellar mass and size 
are of similar strengths. We also demonstrate that there 
is a threshold in stellar mass and size acting at all red- 
shifts which separates galaxies with high crg/S'0.5 and low 
Vrot I 'S'o.s from those with low Ugj 6*0.5 and high Vrot I 'S'0.5 • 

First we consider stellar mass. In Figure |8j as in Fig- 
ure pj we show ag/So,^ and Vrot/ So, ^ as a function of 
redsnift for the mass-limited sample. However, in Fig- 
ure [8] galaxies are divided into three stellar mass bins: 
9.8 < log < 10.1, 10.1 < log M*(M0) < 10.4, 
and 10.4 < log M^Mq) < 10.7. Similar trends to those 
found for the entire mass-limited sample are also found 
separately for the three mass bins. Interestingly, it is 
also evident that at all redshifts examined higher mass 
galaxies are the most kinematically settled (i.e., higher 
Vrot/S^.z and lower (Jg/So.5), intermediate mass galax- 
ies the next most settled, and low mass galaxies the 
least settled. This is consistent with downsizing trends 
in which more massive galaxies att ain final structur e 
and star-formation rates sooner (e.g., |Cowie et al.|1996[ ), 
and we refer to this phenomenon as "kinematic down- 
sizing." Furthermore, there appears to be a threshold at 
log — 10.4 Mq which separates galaxies with high 
o'g/So,^ and low Vrot/So.5 from those with log Cg/S'o.s 
and high Vrot /So. 5- This is most evident in the median 
points in Figure [8] those for the lower and intermedi- 
ate mass bins are very similar, while those for the higher 
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10.4 < log M. (Mq) < 10.7 

10.1 < log M. (UJ < 10.4 
9.8 < log M. (Mg,) < 10.1 




Fig. 8. — As in Figure^ the quantities o'g/S'o.s and Vrot/So.5 are shown as a function of redshift for the mass-Umited sample, but here 
galaxies are divided into three mass bins: 9.8 < log M»(M0) < 10.1 (black squares), 10.1 < log M,(M0) < 10.4 (grey triangles), and 
10.4 < log M«(M0) < 10.7 (red circles). Individual galaxy points are small and medians points in redshift bins are enlarged. Medians 
are shown only for bins in redshift which have > 5 galaxies. AH mass ranges show similar trends with redshift. However, the higher-mass 
galaxies are consistently the most well-ordered at all redshifts (high Vjot/S'o.s and low ag/So.s), the intermediate mass galaxies the next 
well-ordered, and the lowest mass galaxies the least well-ordered. There appears to be a threshold at log M, = 10.4 Mq separating galaxies 
with more ordered motions from those with more disordered motions. 
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Fig. 9. — This figure is the same as Figure Isl except galaxies from the mass-limited sample are divided into three size bins: < 4 kpc 
(black squares), 4—6 kpc (grey circles), and >d kpc (red triangles). Two size measures are considered: Hubble half-light radius which 
is not affected by seeing (top two plots) and emission line extent which is affected by seeing and not corrected for its infiuence (bottom 
two plots). Galaxies in size bins show similar trends to the full mass-limited sample. However, at all redshifts the largest galaxies are the 
most well-ordered (high Vrot/So.s, low (Jg/So.s), the intermediate-sized galaxies the next well-ordered, and the smallest galaxies the least 
well-ordered. There appears to be a threshold at 4 kpc which corresponds to the threshold in stellar mass seen in FigurelT" 
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7. SETTLING FRACTION 

To better quantify how blue galaxies evolve in a-g and 
Vrot^ we use these parameters to define a kinematically 
settled disk galaxy and then study the fraction of settled 
disk galaxies as a function of stellar mass and redshift. 
We return again to the full K07 galaxy sample. For the 
definition of a settled disk, we first look to the distri- 
bution of Vrotl'^g in Figure [3] and then to Hubble mor- 
phologies. Since we regard higher-mass blue galaxies to 
be on average kinematically settled disks at low redshift, 
we look to their values of Vrot/cTg in Figure |3] A typical 
value for galaxies with log M^{Mq) > 9.8 is ^ 3 (i.e., 
log Vrot/cTg = 0.48). We now consider Hubble morpholo- 
gies. K07 and S. Kassin et al. (in preparation) show that 
for the full K07 sample, Hubble morphologies are corre- 
lated with Vrot/o'g such that galaxies with normal disk- 
like morphologies via eyeball inspection have higher val- 
ues than galaxies which appear more disturbed. We find 
that Vrot/f^g ~ 3 indeed provides a reasonable division 
between normal disk-like and disturbed morphologies, as 
demonstrated by Hubble images of 6 galaxies from the 
K07 sample in Figure 



(IFS). 

To compare our dataet with local galaxies, a large 
homogenou s sample is needed, and the GHASP survey 
(Epinat et al. 2010) provides this. Figure 12 shows me- 



dian values of So,5,Vrot, and dg for galaxies in GHASP 
as large open black triangles at 2; = 0. Error bars on 
these points depict the rms scatter of the sample. Al- 
though GHASP spans a larger range in stellar mass than 
our mass- limited sa mple: 9.0 < log M» (Mo) < 11.7 for a 
'diet' Salpeter IMF ( |Epinat e t al.'2008) (which is roughly 
0.15 dex higher than the |Chabri er 2003 IMF used in this 
paper) , the sample medians are consistent with extrapo- 
lations of our relations to z = 0. Individual stellar mass 
measurements for galaxies in GHASP are unavailable to 
create a more comparable sample. 

Kinematic measurements for 47 out of 68 galaxies in 
the IMAGES IFS Survey which have stellar mass mea- 
surements in the same range as our mass- lim ited sample 
are shown as small op en circles in Figure 12 (Puech pri 
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In addition, Vrot/o'g ~ 3 pro- 
vides a good separation~5etwecn the "rotating disks" and 
"perturbed rotators " in the IMAGES Survey at z 0.6 
(|Puech et al.||2007| Figure 6). The median Vrot/cTo for 



2008 
large 



vate communicati on, [ Puech et al. 2008 Neichel et al. 



I Yang et al. 2008). Sample medians are shown as 



lied circles at tiie median redshift of the survey. 



galaxies m the local GHASP Survey (Epinat et al. 12010 1 



is 6.3 with an rms scatter of 3.3, indicating that disk 
galaxies in the local universe are even more kinemati- 
cally quiet than those which pass our "settled" criterion 
This quantitative definition of a settled disk galaxy 



z = 0.6. The sizes used are continuum half-light radii. 
Although IMAGES has a similar selection to our sur- 
vey, the median Vrot and CTg of IMAGES are larger than 
ours at z = 0.6 by 0.23 and 0.18 dex, respectively. The 
measurement methods are sufficiently different for IM- 
AGES compared to our methods that it is difficult to 
compare the two datasets. However, it is possible that 
the IMAGES median cr„ is higher because, unlike our 



(i.e., Vrotlag > 3) can be used to measure the frac- procedure of fitting for K * and a simuhaneously, they 



tion of blue galaxies which are settled, f settle- The evo- 
lution of /settle with redshift for the full K07 sample 
divided into in four stellar mass bins is shown in Fig- 
ure |11[ All bins show increasing /settle with decreasing 
redshift (i.e., increasing /settle with time). For example, 
/settle increases over 0.2 < z < 1.2 from 40% to 89% 
for 10.3 < log M,(M0) < 10.7 and from 31% to 65% for 
9.8 < log M*(Mq) < 10.3. For the less massive bins, we 
probe smaller ranges in redshift. These are lower limits 
to the intrinsic evolution, as discussed in §5.2. Further- 
more, at all redshifts the higher the stellar mass of the 
galaxy population, the higher its value of /settle- The 
same qualitative behavior is found for settled galaxies 
defined as Vrot/<^g > 1 — 4, except for the lowest mass 
bin. 

8. KINEMATIC MEASUREMENTS IN THE LITERATURE 

In this section we compare our mass-limited sample 
with measurements in the literature for local galaxies and 
for galaxies at similar redshifts. We also compile kine- 
matic measurements from the literature and for galaxies 
at higher redshifts. Measurements of 6*0.5, Kotj and Gg 
from the literature are shown in Figure [12] and are de- 
tailed below in order of incre asin g redshift. In the panels 
on the left and right in Figurc[T2l gal axies are color-coded 
according to stellar mass (for a |rhabrier|2003| initial mass 
function) and emission-line extent, respectively, unless 
noted below. The mass-limited sample is also shown in 



are fit for separately in IMAGES. This may result in ag 
measurements which are biased high because pixels with 
the highest signal-to-no ise will be ig,nored when they are 
fit for separately (e.g., ,Davies et al.||2 011[ ). 

The remaining points m l''igiire [12[ plot data for galax- 
ies beyond z = 1.2. Kinematic data are difficult to ob- 
tain at these redshifts since emission lines are redshifted 
to near-infrared wavelengths where the sky is bright in 
continuum and line emission, making spectroscopy diffi- 
cult. Pioneering studies of galaxy kinematics at z > 1 
with single-slit near-infrared spectrographs gave the first 
glimpse of the kinematic state of galaxies at these earlier 
times , albeit only the most highly star-forming systems 



(e.g., Erb et al. 2004 20061. More recently, there has 



been a boom in ll'S observations of galaxy kinematics 
at these redshifts. However, to obtain enough signal-to- 
noise to measure kinematics in a reasonable amount of 
observing time, all samples at these redshifts are unfor- 
tunately still biased toward galaxies with the brightest 
emission lines, which are the most highly star-forming 

systems. 

The first set of points beyond z ~ 1.2 in Figure [12] 
shows kinematic measurements for 29 out of 46 galax- 
ies from the MA S SIV Survey (sn i all op en squares 



tini et al. 2012 Vergani et al. 2012 1. These gala5? 



Con- 



Figure 12 as small grey open circles (labeled "DEEP2" 



with sample medians shown as large filled triangles and 
linear fits to the medians as solid black lines. All liter- 
ature measurements are from integral field spectroscopy 



ies have stellar masses in the same range as our mass 
limited sample and a median redshift of 1.2. Medians in 
'S'o.s, Vrot, and ag for the MASSIV sample are shown as 
large filled squares. The lower mass galaxies in MAS- 
SIV (9.8 < logM*(Mo) < 10.3) have median values 
of these quantities which are consistent with the lower 
mass galaxies in our mass-limited sample. However, the 
higher mass galaxies (10.3 < log M«(M0) < 10.7) have 
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Fig. 10. — Hubble/ ACS V + I images of settled (top) and unsettled (bottom) galaxies which are defined as having Vrot/o-g > 3 and 
Vrot/cTg < 3, respectively. The settled galaxies have normal disk-like morphologies and the unsettled galaxies have disturbed morphologies, 
as determined via eyeball inspection (K07; Kassin et al. in preparation). All galaxies have stellar masses in the range 10.2 < log M« (Mq) < 
11.1, and the images are 6" on a side. 
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Fig. 11. — The fraction of galaxies in the full K07 sample which 
are settled (defined as Vrot/ca > 3), f settle) is shown as a function 
of redshift for four bins in stellar mass. The quantity f settle de- 
creases with redshift (increases with time) for all mass bins. At all 
redshifts the higher the stellar mass of the population, the higher 

its f settle- 

median values which are greater than the higher mass 
galaxies in our mass-limited sample. This is especially 
the case for Kot which has a median value which is 0.20 
dex larger than that for the higher mass galaxies in our 
mass-limited sample at z = 1. This difference is most 
likely due to sample selection, since galaxies in MASSIV 
are chos en to be the very br ightest galaxies in the VVDS 
Survey ( [Contini et aL]|2012[ ). 

Moving to higher redshift in Figure [T2| we note that 
many IFS observations of galaxies at z^l.b have kine- 
matics which cannot be fit by a rotating disk model. 
This is due to a few factors: a non-ordered velocity 
field, small spatial extent, and/or low signal-to-noise. For 
these galaxies we adopt half of the velocity shear across 
the face of the system x Vshear) as Vrot, where Vshear 
is the difference between the minimum and maximum 
rotation velocities and is not corrected for inclination. 
These galaxies are s how n as small filled symbols without 
error bars in Figure [12] to differentiate them from galax- 
ies which were modeled as rotating disks and are shown 
as large open symbols with error bars. 

For 13 extremely high star-for ming galaxies at z ^ 1.3 , 



we adopt measurements from Wisnioski et al. (20111 
which are not modeled as rotating disks (small hlled tri- 
angles in Figure 12). Next, for 8 galaxies at z ^ 1.5, 



we adopt measurements from S. Wright (private com- 
munication) for observations presented in Wright et al. 
( 2007 2009 1 . Galaxies which consist of two or more spa- 
tially separated components are treated as individual sys- 
tems. Four of these galaxies were modeled as rotating 
disks (large open squares) and four were not (small filled 
squares). For galaxies which were so modeled, values of 
Vrot are not corrected for galaxy inclinations. For the 
Wright et al.| ( |2007[ |2009[ ) sample, all but the highest 
value of CTg can be explained by a residual seeing halo, 
and are therefore not likely physical. For an additional 



8 galaxie s at z ~ 1.5, we adopt data from Epinat et| 
al. (2009) and Lemoine-BusseroUe fc Lamareille ( |2010[ ) 
which analyze the same observations. We abbreviate 
these references as "L-BL/E." For both studies, 8 galax- 
ies are modeled as rotating disks. Values of Vrot and ag 
are consistent between the two studies to within uncer- 
tainties except for one gal axy, VVDS 02014710 6. For this 



galaxy we adopt Vrot from Epinat et al. 
model does not extrapolate signihcanti; 
served velocity field. For the remaining 7 galaxies , fitted 



( |2009[ ) since their 
y beyond the ob- 



quantities from Lemoine-Busserolle & Lamareille ( 2010| 
are adopted. Ihese 8 galaxies are shown in I' lgure |12j 
as large open circles. We do not show data for two acT 
dition al galaxies from Lemoine-Busserolle & Lamareille 



(20101 wh ich were not able to be modeled, one of which 
is also in 



Epinat et al 
VVDS 020461235 



(2009). One of these galaxies, 
does not have a measurement of ag . 
The other galaxy, VVDS 020116027, has a velocity field 
which differs significantly between the two studies, likely 
indicating low signal-to-noise. Finally, we show best-fit 
paramete rs for one galaxy at z = 1.6 from IBournaud e"t 



al. ( 2008 ) as large open circles in Figure 12 

For galaxies at z ^ 2, data are adopted from the 
SINS Survey via private communication from N. Forster- 
Schreiber for the most recent models and data as of Au- 

al 



(2008), 



Forster- 
A total 



gust 2011 for g alaxies in [Genzel et 

ihreiber et al.| ( |2009[ ), and |Cresci eTal.i ( |2009| . 
of 14 SUN S galaxies were modeled as rotating disks (large 
open stars in Figure 12 ), and 33 were not (small filled cir 



cles). For galaxies from SINS which were not modeled, 
we estimate ct q as the line-width less the ve locity shear 
(i.e., Vobs from Forster-Schreiber et al.|[2009 less Vshear)- 
For an add i tional 16 galaxies at z ^ 2 we adopt data from 
Law et al. (2009). None of these galaxies were modeled 
as rotating dis ks, and they are shown a s small filled trian 



gles in Figure 12 Three galaxies from Law et al. (20091 
consist of two spatially separated components; we con- 
sider each of these components as a separate galaxy. In 
addit ion, best-fit parameters for on e galaxy at z = 2.03 



from van Starkenburg et al 
open circles labeled "van S." 
For 4 galaxies at z 



(2008) are shown as large 



3, data are adopted from 



Lemoine-Busserolle et al. (2010). Three of these galaxies 
are modeled as rotating disks (large open circles in Fig- 



12 labeled "L-B 



and one is not (filled circle shown 
enlarged in Figure [12] because it overlaps an error bar 
for a different galaxy in the 5o.5 and ag plots). In ad- 
dition, measurements for 19 galaxies at z ~ 3 data are 
adopted from the AMAZE/LSD Survey via private com- 
mu nicatio n from A. Gnerucci for data from IGnerucci eti 
al. (2011 1 and Troncoso et al. (in preparation). A total of 
T2 of these galaxies were modeled as rotating disks (large 
open squares) and 7 were not (small filled squares). 

We will now attempt to summarize the disparate 
datasets from the literature over 1.2 < z < 2.5. Over- 
all, there is large scatter in ag, Vrot, and S'0.5, and no 
clear trends with redshift. While some of the scatter is 
likely due to the range in stellar masses probed, that 
there are no trends is probably also due to other fac- 
tors. Namely, only the brightest galaxies are observed 
at these redshifts, and instrumental effects and poor ob- 
servational conditions likely play roles. It is also clear 
that all high redshift galaxies observed have significant 
contributions to their kinematics from an- 
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Fig. 12. — The quantities plotted here are the same those shown in FigurejS] except here the z axes are expanded to include data from the 
literature. In the plots on the left and right galaxies are color-coded according to stellar mass and emission line extent, respectively, unless 
noted in §8. For the local GHASP survey, only sample medians are shown (large open black triangles) with error bars depicting the rms 
scatter. Large filled symbols denote medians of large samples at z < 1.2: the mass-limited sample from this paper ("DEEP2;" triangles), 
IMAGES (circles), and MASSIV (squares). Individual galaxies from these large samples are shown as small open symbols: DEEP2 (grey 
circles; not colored according to mass/size), IMAGES (circles), and MASSIV (squares). At z > 1.2, only individual galaxies are shown: 
large open symbols with error bars denote galaxies where rotating disk models were fit, and small filled sy mbols without error ba rs denote 
galaxies where models were n ot fit (for these galaxies, Vghear S'ld a-n average ag ar e plotted). Gal axies from'Wisnioski ct al.'2011 ar e shown 
as sn iall filled triangles, from [Wright et al.| 2007 2009 1 as large open squares, from |Lemoine-Bu sscrolle & LamareiUe 2UiO and Epin at~et al.| 
|2009| ( "L-BL/E" ) as large open circles, from Bournau d et al.|2008 as a larg e open circle, from van Starkenburg ct al. 2008 as a lar ge open 
circle, from SINS as large open stars and small filled circles, trom Law et al. 2009 as small filled triangles, trom Lemoine-Busseroll e et al.| 
| |2010[ l ("L-BL") as large open circles and a single filled circle, and trom AMAZE/LSD as large open squares and small filled squares. 
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On average, distant galaxies beyond the last redshift of 
our sample divide into two popula tions (e.g., Law et al. 
2009||Forster-Sclireiber et al |2009[ ): (1) The first consists 
of large massive galaxies with Vrot equal to or slightly 
greater than the largest values which we measure for 
DEEP 2 galaxies, and which have significant values of ct^. 
(2) The second population consists of small lower mass 
galaxies with low or negligible Vrot but with significant 
(T„. These small galaxies are prominent in the Law et al. 



(2009) sample and the SINS survey, but are also found 
aOo^er redshifts in the DEEP2, IMAGES and MASSIV 
surveys. Some of the observations at z > 1.2 may be 
too shallow to detect the full rotation gradients, so until 
these galaxies can be observed with significantly longer 
exposure times to detect possible faint disks, whether or 
not all of them lack significant Vrot is unclear. Interest- 
ingly, at z ~ 3, nearly all galaxies are small but have 
large values of Vrot and Og. 

9. CONCLUSIONS 

We study the internal kinematics of 544 blue galax- 
ies over the last ^ 8 billion years, and find significant 
evolution. Blue galaxies become progressively more well- 
ordered with time as disordered motions decrease and ro- 
tation velocities increase. In addition, galaxy potential 
well depths continuously increase with time. At all red- 
shifts the most massive galaxies are on average the most 
kinematically settled, and the least massive galaxies the 
least kinematically settled. The lowest mass galaxies in 
our survey, which are observed only at low redshift, are 
not well-ordered even today. They may be in the process 
of settling, or may never become settled disks. There 
appears to be a threshold at log — 10.4 Mq and 4.0 
kpc which separates galaxies with more ordered kinemat- 
ics from those with more disordered kinematics. All in 
all, these trends with mass are consistent with downsiz- 
ing trends for other galaxy properties, and we refer to 
them as "kinematic downsizing." 

We define a kinematically settled disk galaxy as having 
a ratio of ordered to disordered motions (K-ot/cg) greater 
than 3. We find that the fraction of settled disk galaxies, 
f settle! increases with time since z — 1.2 for all galaxies 
over 8.0 < log M,(M0) < 10.7. The fraction f settle is 
larger for more massive galaxy populations, independent 
of redshift. These qualitative findings do not change if a 
settled disk galaxy is defined with K-ot/fg > 1 — 4. 

We speculate that the steady kinematic settling seen 
in our data is due to a combination of factors: (1.) The 
frequency of mergers may be decreasing with time. Merg- 
ing, major but more frequently minor, as is expected in 
a ACDM U niverse, will stir galaxies up (e.g 
et al.||2010|). (2. 



Covington 

The amount of mass accreted onto 

galaxies may be decreasing with time. Mass accretion, 
although smoother tha n mergers, might still d i sturb pre- 



existing disks (see e.g ., Bournaud et al. 2011 Cacciato 



Dekel, k, Genel| 2012 for a simple analytic treatment). 



(3.) Galaxies likel y have larger molecular gas reservoir s 
at higher redshift ( IPaddi et al.|2010||Tacconi et al.|20To| . 
Higher gas fractions should result in more star-formation, 
and feedback from star-formation may a lso stir up the 
gas in galaxies (e.g.. Silk fc Norman|2009 and references 
therein). (4.) Finally, higher gas tractions can also lead 
to violent disk instabilities 



which increase random motions, and may also be a step 
in the direction of more star-formation. 

At bottom, the primary causative factors of disk set- 
tling are likely two-fold: more merging/accretion and 
higher gas fractions at early times. Faster star-formation 
rates are likely byproducts of these two factors, and 
increasing potential well depths a byproduct of merg- 
ing/accretion. Since both factors decline with time, a 
general kinematic settling with time is expected. Fur- 
thermore, both factors are apparently declining earlier 
in massive galaxies, at least since 2^1, providing a 
natural explanation of the kinematic downsizing we find. 

It is yet unknown whether disk galaxies in cosmolog- 
ical numerical simulations undergo a kinematic settling 
akin to what we find in this paper. Only recently have 
hydrodynamic simulations of galaxy formation been run 
with enough numerical resolution to study the effects of 
mass accretion in detail. Simulations find that mass ac- 
cretes directly on t o galaxies from cosmic filaments (e.g., 
Katz et al.|[2003l [Birnbo i m fc Dekel||20031 |K eres et al. 
k et al 



2005| lOcvirk et al.||2008[ [Brooks et al.|l2009[ ). The ef- 
fects of accretion have been studied wi th small samples 



of si mulated galaxies at z ~ 2 (e.g., [Bournaud et al. 
2669;' Ceverino et al.|2010 ). These studies find that mass 
accretion perturbs disks and leads to violent disk insta- 
bilities, but the effects of accretion have yet to be studied 
for a representative sample of simulated galaxies at lower 
redshifts. The simulations which have been run to lower 
redshifts generally focus on reproducing th e photometric 
and structural properties of galaxies (e.g., Governato et 



|al.l2007[[2009[ [Brook et al.|2011||MaFtig et al.|^Ol^[ and 

references therein), and not (Ja, Vrot, an d iS'0.5 . On the 



Cacciato, Dekel, & Genel 2012 and references therein) 



Bournaud et al. 2011 



analytic side, Cacciato, Dekel, & Genel (2012) created 
a simple model tor the kinematic evolution 01 massive 
disk galaxies at z ~ 2. They predict that these galaxies, 
if still around today, are massive disks with similar ro- 
tation velocities but less disordered motions, similar to 
the kinematic settling we find. Numerical simulations are 
not yet ripe for studies of the kinematic evolution of blue 
galaxies since z ~ 1 , and we look forward to comparisons 
with future predictions. 

Our data are augmented with a compilation of mea- 
surements of the kinematics of emission-line galaxies 
from the literature over < z < 3.8. Observations of 
local massive disk galaxies reveal few examples where 
disordered motions (as quantified by Ug) play as an im- 
portant role as in si milar-mass eniission -line galaxies at 
higher redshifts (e.g., Epinat et al.|2010 ). Measurements 
of galaxy kinematics at z > 1.5 show few clear trends 
with redshift, likely primarily due to sample selection 
since only the most highly star-forming systems at these 
redshifts can currently be studied, but also due to instru- 
mental artifacts and observational conditions. However, 
it is clear that higher redshift galaxies have a significant 
amount of disordered motions, whether or not they show 
evidence of rotation. 

In summary, galaxies seem to have a life cycle: Early in 
their lives they are accreting baryons rapidly, undergoing 
mergers, and possess a large amount of gas. Later in their 
lives the accretion decreases along with their gas content. 
Kinematic settling appears to be a natural, additional 
facet of this basic evolutionary arc. 
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